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Abstract A particle-in-cell code is used to examine 
contributions of the pickup ions (Pis) and the so- 
lar wind ions (SWs) to the cross shock electric field 
at the supercritical, perpendicular shocks. The code 
treats the pickup ions self-consistently as a third com- 
ponent. Herein, two different runs with relative pickup 
ion density of 25% and 55% are presented in this pa- 
per. Present preliminary results show that: (1) in the 
low percentage (25%) pickup ion case, the shock front 
is nonstationary. During the evolution of this perpen- 
dicular shock, a nonstationary foot resulting from the 
reflected solar wind ions is formed in front of the old 
ramp, and its amplitude becomes larger and larger. At 
last, the nonstationary foot grows up into a new ramp 
and exceeds the old one. Such a nonstationary pro- 
cess can be formed periodically. When the new ramp 
begins to be formed in front of the old ramp, the Hall 
term mainly contributed by the solar wind ions becomes 
more and more important. The electric field E x is dom- 
inated by the Hall term when the new ramp exceeds 
the old one. Furthermore, an extended and station- 
ary foot in pickup ion gyro-scale is located upstream 
of the nonstationary /self-reforming region within the 
shock front, and is always dominated by the Lorentz 
term contributed by the pickup ions; (2) in the high 
percentage (55%) pickup ion case, the amplitude of the 
stationary foot is increased as expected. One striking 
point is that the nonstationary region of the shock front 
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evidenced by the self-reformation disappears. Instead, 
a stationary extended foot dominated by Lorentz term 
contributed by the pickup ions, and a stationary ramp 
dominated by Hall term contributed by the solar wind 
ions are clearly evidenced. The significance of the cross 
electric field on ion dynamics is also discussed. 

Keywords Pickup ions; nonstationary shock; Particle- 
in-cell simulation 



1 Introduction 

Collisionless shocks are of great interests in space 
physics, plasma physics and astrophysics since many 
theoretical, numerical and experimental studies have 
already evidenced that in the shock transition the bulk 
flow energy of p lasmas is conv e rted irreversibly into 



thermal energy (ISagdeevI Il966l ; ILembege et al.1 12004 ; 



Burgess et all 120051) . During the conversion process, 
the interaction between the electromagnetic fields and 
particles replaces the role played by collisions in a 
normal hydrodynamics. In quasi-perpendicular shocks 
{Obu > 45°, 9bu is the angle between the upstream 
magnetic field and shock normal), the normal/cross 
shock electri c field / potential plays an key role in both 
acceleration ( Zank et al.lll996l ; IShapiro and Ucei 2003 ; 



Yang et al. 2009a) and thermalization (IBurgess et al 



Il989t ILembege and Savoinilll992l) of the incident parti- 
cles. In-situ observations of the electric field at Earth's 
bow shock provide important and unique experimen- 
tal data that can be used to understand the details 
of plasma thermalization and energetic particle pro 



ductio n in more general circumstances. IWalker et al 
( 20041 ) experimentally studied the structure of the elec- 
tric field profile during a number of crossings of the 
quasi-perpendicular bow shock using data from the 
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Cluster satellites. Their results showed that the am- 
plitude of the cross shock electric field appears to in- 
crease as 0Bn approaching 90°. Based on Cluster four- 
spacecraft observations, contributions to the electric 
field of the terrestrial bow sh ock are analyzed in de- 
tail by lEastwood et all (|2007t ). They found that the 
Hall term in the generalized Ohm's law accounts for 
a majority of the normal electric field of supercritical, 
quasi-perpendicular shocks by two different methods: 

1. the single-spacecraft technique in which the shock is 
assumed as one-dimensional (TP) and stationa ry, and 

2. the curlometer technique ([Dunlop et al.ll2002T ) which 
requires that the spacecraft are co-located in the shock 
ramp and the shock is stationary among the satellites. 

Both particle-in-cell(PIC) and hybrid simulations 
have clearly evidenced that supercritical, quasiper- 
pcndicular shocks are nonstationary and suffer a self- 
reformation on the loca l gyro s cale of the incide n t ions 



( Lembege and Dawsonl Il987j: Nishimura et al 



2003; 



field is dominated by the off-diago nal electron pressure 
|Vasvliunaslll976t IPritchettl lioQ]}) . In order to under- 
stand the physics in different scales (e.g. the whole 
shock front is in ion inertial scale and the ramp is in 
electron inertial scale) of the shock front struct u res of a 
supercritical, perpendicular shock, I Yang et al.l (|2009bl ) 
borrowed the same techniques which has been success- 
fully applied to the research on collisionless magnetic re- 
connection. First, a PIC simulation including all of the 
behaviors of ions and electrons were performed. Non- 
stationary fields and corresponding particle dynamics 
were obtained self-consistently. Second, contributions 
of Lorentz, Hall and electron pressure terms to the 
cross shock electric field were calculated by using the 
self-consistent particle data issued from the PIC sim- 
ulation. The main conclusions are the following: (1) 
the contributions of Lorentz, Hall and electron pressure 
terms to the cross shock electric field is evolving ver- 
sus tim e. The electron pressure term within the shock 



Scholer and Matsukiyo 2004; Hcllinge r et al. 2007; Lembegcfntnihjways contributes little; (2) when the new ramp 
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^ Yang et all2011a ; Shinohara et al.ll2011 ). Hada et al-^ egins to be formed in front of the old ramp, the Hall 



(2003) analyzed the time evolution of the normal elec 
trie field and the fraction of reflected ions at the re- 
forming shock front by using both PIC simulation and 
semi-numerical model. Their findings suggest that the 
reflected ions have a very coherent motion (narrow ring) 
as described by the so called trapping loop/vortex in 
the phase space (V x — X) at the beginning of the reflec- 
tion. This coherent forces these ions to accumulate over 
a narrow spatial range at a foot length distance from 
the ramp; it results that ramp and foot are well sepa- 
rated. This accumulation is suitable for building up a 
peaked-like normal electric field which itself reflects at 
later times new incoming ions at the foot before these 
reach the ramp. The peak of electric field produces an 
enhancement of local magnetic field. At last, a new 
ramp (i.e. grown-up foot) is fully built up and sub- 
stitute the old one. This process is periodic, and is 
named self-reformation. Nonstationary characteristic 
of the shock front is later confirmed by an individual 
case st udy of Cluster obse r vation at the t e rrestr ial bow 



term becomes more and more important; (3) the elec- 
tric field E x is dominated by the Hall term when the 
new ramp exceeds the old one. Nevertheless, the shock 
front topology of a nonstationary shock can be modified 
by the local pickup ions (Pis). 

Pickup ions are extensively and directly observed 
in the solar wind upstream of the Earth's bow shock 
(at ~1 AU) wit h the SULEICA instrument on aboar d 
AMPTE/IRM jMobius et al.l Il985t lOka et al.l l2002h . 
around the shock associated with the corotating in- 
teraction region (CIR) in the interplanetary space (at 
~4.5AU) w ith SWICS instrumen t onboard the Ulysses 
spacecraft (jGloeckler et al.lll994l) , and in the upstream 
of the Halley comet bow shock with HERS instr ument 
on the Giotto spacecraft ( Neugebauer et alll 19891) . The 
production process and feature of the pic kup ions in the 



solar wind have been s tudied quite well (Sa gdeev et al 



1986; iLee and Iplll987l ). Neutral atoms and molecules, 



shock (jLobzin et al.ll2007t ). iMazelle et "all fj2010h statis- 



tically confirms this robust characteristic. In addition, 
they found the magnetic field ramp of the shock often 
reaches a few c/w pe , and the time variability of scales 
observed for the same upstream conditions is also a 
clear signature of the nonstationary shock front. 

In recent study on collisionless magnetic reconnec- 
tion. First, at the ion inertial length scale, ions an d 
electrons are decoupled ( Pritchettll2001 ; Lu et alj2011 ). 
The resulting Hall effect determines the reconnection 
electric field. Second, at the electron inertial scale, 
even the electrons are demagnetized, and the electric 



penetrating from the local interstellar medium or escap- 
ing from comets, are ionized by photoionization, elec- 
tron impact or charge exchange with the solar wind. 
Under the influence of the solar wind magnetic field 
the ions compose a ring-beam distribution in velocity 
space. Then they are scattered by ambient and excited 
Alfvenic fluctuations to form a spherical shell distribu- 
tion. Interstellar neutral hydrogen (H + ) has a much 
higher density than other interstellar neutrals; inter- 
stellar helium accounts for ~ 10% of the interstellar 
gas and its ioniza tion rate is ~ 10% that of hydro- 
gen riMobiusl ll986). The percentage of pickup ions to 
the solar wind ions (SWs) at 50 AU is estimated to be 
10% and increase linearly with distance from the sun 
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( Vas vliunasl 1 1 9 761 ) . So the pickup H + is expected to be 
the dominant (i.e. number density of pickup protons is 
> 20% that of the solar wind protons) pickup species 
at the termination shock (TS) at abou t 70 ~ 100 AU 
(|Richardson et al lboodlWu et al.ll2009l). Early 1-D hv 



Table 1 Upstream plasma parameters defined for the 1-D 
PIC simulation 



brid simulations ( Omidi et al.lll986t lOmidi and Winske 



1987t iLiewer et al.l 119931 ) have revealed that an ex- 
tended foot, which scales in reflection ion gyroradius, is 
seen on the supercri tical shocks as a results of th e re- 
flected pickup ions. iMatsukivo and Scholerl (|201ll) ex- 



amined the effect of pickup ions on the heliospheric ter- 
mination shock by using a 1-D PIC code. They found 
the pickup ions in the extended foot causes an electric 
field in the shock normal direction. This leads to a large 
increase of the shock potential barrier upstream of the 
magnetic field ramp. However, the impact of pickup 
ions on the self-reformation process and contributions 
to the cross shock electric field at nonstationary shocks 
has not been analyzed yet. This gave us the motivation 
for investigating the impact of Pis on the relative con- 
tributions of Pis and solar wind ions (SWs) respectively 
to the total normal electric field of the shock front. 

Th is paper represent s an extension of a previous 
work ( Yang et al.l2009bl) . Herein, we aim to address the 
following questions by using one-dimensional PIC simu- 
lation (the code treats the pickup ions self-consistently 
as a third component): (1) Which term (such as Hall 
term and Lorentz term) in the generalized Ohm's law 
accounts for a majority of the normal electric field of 
a supercritical, perpendicular shock in the presence of 
pickup ions? Does it vary versus time? (2) What is 
the contributions of pickup ions and solar wind ions to 
the terms mentioned above? (3) How is this impact ac- 
cording to different percentages of upstream pickup ions 
on this competition? (4) What are the consequences of 
high percentage of pickup ions on the nonstationarity of 
the shock front? This paper is organized as follows. In 
section 2, we briefly describe the numerical model used 
herein. Sections 3 and 4 present the simulation results 
at (a) low Pis percentage case, and (b) high Pis per- 
centage case, respectively. The main conclusions will 
be summarized in section 5. 



2 Simulation Conditions 

We use a one-dimensional electromagnetic PIC code 
to simulate the evolving structure of a supercriti- 
cal, collisionless, perpendicular shock. Simulations 
of nonstation ary shocks have b een performed, for 
example, by lYang et al.l (|2011bf ). where the pickup 
ions (H + ) arc treated as test particles. Herein, 
the shock is produced by the injection method as 
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in previous wor ks dQuesd Il985l burgess et al.l Il989i: 
Nishimura et al.l 120031: IChapman et al.l l2005h . In 



or- 
der to self-consistently reproduce the impacts of pickup 
ions, we simulate three species within our code: elec- 
trons, solar wind ions, and pickup ions. Particles are 
injected on the left-hand side of the simulation box with 
a inflow/upstream drift speed Vi„j, and are reflected at 
the other end. The distribution functions for the so- 
lar wind ions and electrons are Maxwellian. Pickup 
ions are injected on a thin sphere in velocity space 



centered at Vj„j with r 


adius V s hM as in earlier works 


(Kucharck and Scholer 


1995: Lee et al. 


2005: 


Wu et al. 



by an ion/ion beam instability (|Scholer and Matsukivo 



2004}) and moves with a speed V re f from the right-hand 
side toward the left. The upstream Alfenic Mach num- 
ber of the shock is Ma — (Vi n j + V re f)/VA , where the 
Alfcn speed Va is equal to 1 . All basic parameters are as 
follows: Plasma box size length L x — 92.16(c/cj p i); ve- 
locity of light c = 21.4, and mass ratio mi/m e — 1836. 
The electron/ion temperature ratio T e /Ti = 1.36 is cho- 
sen in order to investigate the particle acceleration at 
a reforming shock. The ambient magnetic field along 
the y direction is \E>o\ = 1- Low beta values /3, = 0.125 
and /3 e =0.1 are chosen in order to investigate the par- 
ticle acceleration at a reforming shock. Initially there 
are 200 particles of each species in a cell. The up- 
stream plasma is quasi- neutrality, i.e. n e = n;, where 
rii = nsws + npis = N S ws x SWs% + N PIs xPIs%. n e , 
fii, nsws and npj s are densities of the electrons, the to- 
tal ions, the solar wind ions and the pickup ions, respec- 
tively. Nsws, Npi s , SWs% and PIs% are the counts 
of each species of ions and their relatively weighted per- 
centage in the PIC simulation. For these initial condi- 
tions, the plasma parameters are summarized in Table 
Q]for both electrons and SWs. The Larmor gyro radius 
in the table is calculated based on the thermal velocity. 
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3 Simulation Results 

This section can be separated into two parts. In the 
first part, we show the contributions of solar wind ions 
and pickup ions to the cross shock electric fields at three 
typical shock profile issued from the low PIs% case. In 
this case, the percentage of the pickup ions is 25%, and 
the shock front is characterized by a self-reformation 
process with a cyclic period about 1.36fi7^ (f2 c j is the 
upstream ion gyro- frequency) . During this reformation, 
the instantaneous Mach number of the shock at differ- 
ent times can differ from the average value 71^4=5. 29. 
In the second part, we do similar things for high PIs% 
case, in which the percentage of the pickup ions is set 
to 55%. In this case, the shock front nonstationarity 
is fully suppressed. In each part, the analysis is based 
on systematic comparison between the Hall term and 
Lorentz term (contributed by solar wind ions and/or 
pickup ions) in the generalized Ohm's law in order to 
clarify the answers to the questions addressed in this 
paper and mentioned in section 1. 

3.1 Effect of Low Percentage Pickup Ions on the Cross 
Shock Electric Field (nonstationary case) 

Firstly, let us analyze the impact of low percentage 
pickup ions on the cross shock electric filed at a non- 
stationary shock. The nonstationarity of the perpen- 
dicular shock can be demonstrated in figure 1, which 
shows the time evolution of the magnetic field B y from 
t =4.6 to 120 ( ~ 1 . The shock is propagating from the 
right to the left. At about t =4.6, the shock ramp 

("old" ramp) is at about X =84.48c/Q pi . Later at temperature (jrleinz et al.ll2005l) 



Lee et alj 120051: IMatsukivo et al]|2007t IWu et al.ll2009 ; 
iMatsukivo and Scholerl 2011 ). In this subsection, we 
will clarify why the low percentage pickup ions can not 
suppress the self-reformation of the shock front by an- 
alyzing the cross shock electric field which play a key 
role in the reflection of incident ions. 

Let's concentrated on the evolution of the electric 
field in this nonstationary perpendicular shock. Assum- 
ing a two-fluid model for the plasma, the electric field 
can be obtained directly from the electron momentum 
equation without further approximation 



dt 



(1) 



Here U e is the electron fluid velocity, B the magnetic 
field, P e the electron pressure tensor, m e the electron 
mass, and e the electric charge. If we use the total 
current density J = e{iii\Ji — n e \J e ) (where Ui is the 
ion fluid velocity) and considering the one-dimensional 
property of our simulation, the normalized Eq. ([T]) can 
be rewritten as 
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(2) 



where T P , 



J2 \Vex,k 

fe=i 



U ex ) 2 /n e 



^exy Z^/ i^ex,k 
k=l 



U ex ) x {V ey ,k ~ U ey )/n e , and T exz = J2 (K: 
(V e: 



U ex ) X 



fe=l 



Uez ) I n e are components of the electron kinetic 



about t =5. 2D" 1 , the foot amplitude becomes larger 
and larger due to the accumulation of the reflected so- 
lar wind ions in the foot. Then a "new" shock ramp 
is formed. At about t =5. 651,7/, it exceeds the am- 
plitude of the "old" ramp, and the "new" shock ramp 
is around X =81.89c/f2 pi . Simultaneously, the "old" 
shock ramp becomes weaker and weaker, and is lo- 
cated downstream of the "new" ramp. The shock 
front is characterized by a self-reformation with a cyclic 
period about 1.36f2~ 1 . The self-reformation of the 
shock front is due to the coupling between the "in- 
coming" a nd "reflected" ions , which has been demon- 
strated by lHada et al.l (|2003h at supercritical perpen- 
dicular shocks. We see that the solar wind ion dynamics 
proceed essentially as in the absence of the pickup ions, 
whereas the pickup ions form an stationary extended 
foot starting ~ 4(c/f2 C i) upstream. This phenomenon 
has already been demo nstrated by a good deal of hybrid 
and PIC simulations ( Omidi et al.l Il986t iLiewer et al 



19931: iLipatov and Zankl Il999t IChapman et al.l 12005 



x , Gy , and 6^ are 
unit vectors. The first, second, third and fourth terms 
on the right of Eq. ([2]) are Lorentz term, Hall term, 
electron pressure gradients term and electron inertial 
term, respectively. All of the physical quantities on the 
right of Eq. ([2]) can be statistically obtained from the 
PIC simulation. We choose three typical shock pro- 
files at three different times (marked by "A" , "B" and 
"C" in figure 1) within one shock self-reformation cycle 
(from t =10.ir2~ 1 to 11.4ft ) and analyze separately 
the corresponding roles of Lorentz and Hall terms in 
the cross shock electric field E x . Herein, we concen- 
trate on the competition of the solar wind ions and 
pickup ions for Lorentz and Hall terms during the shock 
front self-reformation. The electron pressure gradients 
term and electron inertial term are neglected because 
they do not play impo rtant roles in ion r eflection and 
shock nonstationarity (jYang et al.l l2009bl ). Since the 
normalized bulk velocities of the total ions, SWs and 

Nsws N PIs 

Pis are U 4 = ( £ V S Ws,k X SWs% + £ V PIs , k x 



k=l 



k=l 
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Nsws 

PIs%)/nj, U S ws = E V S Ws,k/nsws and U P/s = 
fc=i 

N PIs 

^2 "V Pis,k/npi s , respectively. The normalized cur- 

fc=i 

rent densities contributed by the total ions, SWs and 
Pis can be written as J = rijUj — n e U e , Jsws = 
nswsUsws - nsWs n e U e and Jp/ S = npi s X3pi s - 



npi s 



e U e , respectively. Eq. ([2]) becomes: 



-Ex = (U S ws x B) x (U P/S x B)^ + 

n e n e 

— (Jsws x B) x + — (Jp/ S x B) x , (3) 

n e n e 

where the first and second terms on the right of Eq. ([3]) 
are contributions of SWs and Pis to the Lorentz term 
of E x , and the other two are contributions of SWs and 
Pis to the Hall term of E x . 

Figure 2a-2c shows the snapshots of CSP (cross 
shock potential), B y and E x (in black) at shock pro- 
file A (at t = 10. 4f2~ 1 ) of a reforming shock, re- 
spectively. At this time, the shock front can be di- 
vided into three regions: region I (highlighted in cyan) 
is the stationary foot dominated by the pickup ions, 
region II (highlighted in yellow) indicates the time- 
varying/nonstationary foot results from the accumula- 
tion of the reflected solar wind ions, region III (high- 
lighted in magenta) is the nonstationary shock ramp. 
The position of the ramp "O" (denoted by vertical 
dashed) is at X —71.85c/uj p i. The width of the shock 
ramp is measured from the maximum point of the 
cross shock electric field —E x to the maximum point 
of the magnetic ov ershoot (jShimada and Hoshino|[2005 : 
Yang et al.ll2011al) . is about 0.l5c/uj p i =6Ac/uj pe . And 
the widths of foots in region I and II are about Ac/w P i 
and 1.3c/wpi, respectively. Figure 2a shows the phase 
space plots {Vi X — Xj) of the solar wind ions. A frac- 
tion of the them are reflected at region III and accumu- 
lated at region II. Corresponding phase space plot of the 
pickup ions is not sho wn here, and it has alrea d y bee n 
analyzed in detail by IMatsukivo and Scholerl ([20111 ). 
On this point, there is no longer to do more. Figure 
2b shows the density profiles of the Pis (blue) , the SWs 
(red) and the total ions (green). In region I, the con- 
tribution of Pis is about two-thirds of that of the SWs. 
The density profile of SWs in this region is flat and 
the density of Pis increases slightly. In region II, the 
density of SWs increases distinctly, and the density of 
Pis becomes flat. And in region III, the local density 
profile of SWs increases precipitously because most of 
the incident solar wind ions are reflected at the ramp, 
and the density profile of the Pis is still flat. In fig- 
ure 2c, the contributions of the Pis, SWs and the total 



ions to the Hall term of the cross shock electric field 
E x (black curve) are indicated by blue, red and green 
curves, respectively. The main features are summarized 
as follows: 1. in region I, the amplitude of the total Hall 
term is close to zero. Because the contributions of Pis 
and SWs to the total Hall term (in green) are roughly 
equal and opposite to each other. 2. in region II, the 
Hall effect of the Pis (i.e. the term ^-(Jp/ S x B) x in 
Eq. ([3]), in blue) is almost zero, and the Hall term is 
dominated by SWs. 3. in region III, the Hall effect 
associated with the SWs is much larger than that of 
the Pis. The spatial form of the Hall term contributed 
by SWs (i.e. the term ^(J SH / S x B) x in Eq. ©, 
in red) apparently matches with that of E x . In figure 
2d, the contributions of the Pis, the SWs and the to- 
tal ions to the Lorentz term of the cross shock electric 
field E x (black curve) are color-coded severally as in 
figure 2c. The main features are as follows: 1. in re- 
gion I, the effect of SWs on the total Lorentz term (in 
green) is very small. The Pis dominated the Lorentz 
term in this region. 2. in region II, the contributions 
of SWs (i.e. the term (U S w s X B) x in Eq. ©, 

in red) is a bit smaller than that of Pis (i.e. the term 
_n£i»(Upj s x B) x in Eq. ©, in blue) to the total 
Lorentz term. 3. in region III (i.e. in the vicinity of 
the shock ramp) , the contribution of Pis to the Lorentz 
term is almost zero. 

In summary, the cross shock electric field E x in re- 
gion I (i.e. at the stationary extended foot) is domi- 
nated by the Lorentz term, and is mainly contributed 
by Pis. In region II (i.e. at the time-varying solar 
wind ion foot, which will be clarified in the remain- 
der of this subsection), the contributions of the Lorentz 
term and Hall term to the total E x are comparable to 
each other. And the profile E x in region III (i.e. at 
the steep ramp) is dominated by the Hall term, and is 
mainly contributed by SWs. 

Figure 3a-3c shows the snapshots of CSP (cross 
shock potential), B y and E x (in black) at shock pro- 
file B (at t = 10. ) of a reforming shock, respec- 
tively. At this time, the whole shock front still can be 
divided into three regions as in figure 2. The amplitude 
of the old ramp "O" is decreased, and its width is about 
0.13c/uj p i =5.6c/w pe . It located at X =70.88c/w P j. 
The extended foot caused by pickup ions in region I is 
stationary, and the nonstationary foot caused by solar 
wind ions in region II becomes broader. Their widths 
are about Ac/uj p i and 1.55c/w P i, respectively. Now, 
the cross electric field E x still has a negative value at 
the shock old ramp ("O"). However, its amplitude is 
strongly decreased and the foot in region II grows into a 
new ramp ( "N" ). Figure 3a shows the phase space plots 
(Vix ~ Xi) of the solar wind ions. It clearly appears that 
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the reflected ions have a very coherent motion (within 
regions II and III) as described by the well defined trap- 
ping loop (vortex). This vortex is a manifestation of 
shock front self-reformation [Hada et al., 2003]. At this 
time, the bulk velocity of the incident solar wind ions 
begin to decrease at the new ("N") ramp in region II 
instead the old one in region III. Figure 3b shows the 
density profiles of the Pis (blue), SWs (red) and the 
total ions (green), as in figure 2b. In region I, the total 
density profile gradually increases versus X due to the 
accumulation of the reflected pickup ions. More and 
more solar wind ions, which are reflected at the ramp 
in early time, are accumulated in region II. Thus a new 
ramp ( "N" ) appears and its amplitude reaches at about 
87.5% of the old ramp ("O"). The reflected solar wind 
ions feeds up the self-reformation of the shock front. 
And in region III, the density profile of the SWs is de- 
creased. The local electromagnetic fields are weaker 
than that at profile A. The density profile of SWs in 
region II and III have a relationship of restricting each 
other. And the density profiles of Pis in region II and 
III are still flat as that in figure 2b. The main features 
shown in figure 3c are summarized as follows: 1. in 
region I, the amplitude of the total Hall term is close 
to zero, as in figure 2c. 2. in region II, the spatial form 
of the Hall term contributed by SWs (red) fine matches 
with that of the total Hall term. 3. in region III, the 
Hall term is still dominated by SWs, but its amplitude 
is decreased. The main features of figure 3d are sum- 
marized as follows: I. in region I, the Pis still dominate 
the Lorentz term as in figure 2d. 2. in region II, the 
contributions of SWs to the Lorentz term is as much as 
that to the Hall term in figure 3c. 3. in region III, the 
contribution of Pis to the Lorentz term is almost zero. 

In summary, the cross shock electric field E x in re- 
gion I (i.e. at the stationary extended foot) is domi- 
nated by the Lorentz term as at profile A, and is mainly 
contributed by Pis. In region II (i.e. at the nonstation- 
ary solar wind ion foot), the Lorentz term and Hall 
term are comparable to each other, and both of them 
are increasing. The amplitudes of the total E x and the 
cross shock potential in region II have already exceed 
that in region III. Furthermore, Hall term in this region 
is dominated by SWs. In region III, the amplitude of 
E x is decreased a lot. 

Figure 4a-4c shows the snapshots of CSP (cross 
shock potential), B y and E x (in black) at shock pro- 
file C (at t = 11.3f2^ ) of a reforming shock, respec- 
tively. In contrast to figure 2 and 3, a striking point 
is that the whole shock front only can be divided into 
two regions because the amplitude of the new ramp 
("N") have already fully overcome the old one at this 
time. The new ramp "N" in region III is strong enough 



to reflect the incoming upstream solar wind ions, and 
its width is about 0.1&c/u pi =6.9c/cj pe . It located at 
X =69.48c/w p i. The extended foot in region I is still 
stationary, and its width is about 4c/cj pi . Different 
from the results shown in figure 3, figure 4c clearly evi- 
dence that there is no region II foot ahead of the steep 
ramp, and the Hall term dominated by SWs perfectly 
matches with the total cross shock electric field E x in 
region III. Figure 4d shows the Lorentz term dominated 
by Pis apparently matches with the E x in region I. 

3.2 Effect of High Percentage Pickup Ions on the 
Cross Shock Electric Field (stationary case) 

Secondly, let us analyze the impact of high percent- 
age pickup ions on the cross shock electric filed at a 
supercritical perpendicular shock. In the current case, 
the percentage of the pickup ions is set to 55%, and the 
other parameters of the simulation are kept unchanged. 
In contrast to the above results (section 3.1), a sta- 
tionary shock can be demonstrated in figure 5. This 
stationary, supercritical shock is characterized by an 
overshoot, a ramp and broad foot in front of the ramp. 
Without loss of generality, we analyze the cross shock 
electric field E x at t ^IQAVt^ 1 (marked by "A"' in fig- 
ure 5) because the shock structure is not changing with 
time. 

Figure 6a-6c shows the snapshots of CSP (cross 
shock potential), B y and E x (in black) at shock pro- 
file A' (at t — 10.417,7/) of a reforming shock, respec- 
tively. At this shock profile, the whole shock front can 
be divided into two parts: region I and region III simi- 
lar to that of shock profile C mentioned in section 3.1. 
The shock ramp is at X —66.58c/u} pi , and its width is 
about 0.l9c/uj p i =8.2c/oj pe . Compare with the results 
described in section 3.1, the amplitudes of the overshoot 
and cross shock electric field at shock profile A' (i.e. at 
a stationary shock) is much lower than that measured 
at a reforming shock, e.g. at profiles A (Figure 2) and 
C (Figure 4). Figure 6a shows the solar wind ions re- 
flected at the shock ramp (marked by vertical dashed 
in region III) can not form a clear vortex which is used 
as a key evidence of shock front self-reformation [Hada 
et al., 2003]. Figure 6b shows the density profiles of 
the Pis (blue), SWs (red) and the total ions (green), 
as in figure 4b. In region I, the total density profile 
is increased, and its shape is mainly controlled by Pis. 
The density of SWs exceeds that of Pis in region III. 
Figure 6c and 6d shows the contributions of SWs and 
Pis to the Hall term and Lorentz term at the station- 
ary shock (profile A') are analogous to that obtained 
at the self- reforming shock (e.g. profile C in figure 4). 
The only difference is that the amplitude of the cross 
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shock electric field E x in region III is much lower, and 
the foot in region I is slightly higher. 

In summary, high percentage Pis leads to a station- 
ary shock (i.e. without self-reformation) even in low 
case. This point al r eady h as been pointed out by 



Matsukivo and Scholerl (|201ll ) in which they claimed 
that a shock becomes stationary even with 30% Pis. 
The cross shock electric field E x in region I is still dom- 
inated by the Lorentz term contributed by Pis. In re- 
gion III, E x is still dominated by the the Hall term 
contributed by SWs. These two points are the same 
as that obtained at the steep and narrow shock profile 
(e.g. profiles A or C) in low percentage Pis case (i.e. 
nonstationary case). Furthermore, low electromagnetic 
field within the shock front can not sustain the ion vor- 
tex in phase space. Thus the region II with feature of 
self-reformation disappears at this stationary shock. 



4 Conclusions 

In this paper, we used 1-D PIC simulations to study the 
contributions of the cross shock electric field in super- 
critical, perpendicular shocks in the presence of pickup 
ions. Firstly, we separate the electric field into several 
terms (such as Hall term and Lorentz term), and their 
importance is evaluated. Secondly, we compare the con- 
tributions of pickup ions and solar wind ions to each 
term above. Finally, we study the impact the percent- 
age of pickup ions on the shock front nonstationarity. 
The analysis has evidenced the following features: 

1. In low PIs% case, the supercritical shock is non- 
stationary. For the electric field E x , after the new shock 
ramp is formed, the Hall term becomes more and more 
important with the increase of the new ramp. After 
the new ramp exceeds the old ramp, its width is very 
small, which is comparable with or even smaller than 
the ion inertial length, and the Hall term dominates the 
electric field E x . At that time, the electric field E x is 
obviously larger that a t other times. As pointed pre- 
viously by lYang et al. ( 2009al ) the shock surfing accel- 
eration (SSA), where the particles are reflected mainly 
due to the electric field E Xl ion Lorentz term become 
negligible at that time. The characteristics mentioned 
above is nearly the same as that obtained at the non- 
stati qnary perpendicula r shock in the absence of pickup 
ions (|Yang et al.ll2009bl) . 

2. In low PIs% case, the extended foot caused by 
the reflected pickup ions also has been retrieved. E x in 
the extended, stationary foot region (so-called region 
I) is mainly contributed by the Lorentz term, which is 
always dominated by the pickup ions. In solar wind 
ion foot region (so-called region II), contributions of 



the pickup ions and solar wind ions are competing with 
each other in the Lorentz term of E x . Simultaneously, 
the Hall term is mainly contributed by the solar wind 
ions. And in the ramp and overshoot (so-called region 
III), E x is mainly contributed by the Hall term and is 
always dominated by the solar wind ions. 

3. In high PIs% case, the amplitude of the extended 
foot is increased as expected. The striking point is that 
region II can not be found any more. The whole shock 
front becomes stationary. Because the amplitude of the 
overshoot and cross shock electric field is decreased and 
thus they are not strong enough to form the solar wind 
ion vortex in phase space at the ramp. 

4. In high PIs% case, there are two remainder re- 
gions (I and III) within the shock front. E x in region 
I is dominated by the Lorentz term contributed by the 
pickup ions and E x in region III is dominated by the 
Hall term contributed by the solar wind ions, respec- 
tively. 

Recently, although the real ion-to-electron mass ra- 
tio has been employed. A few other plasma parameters 
are demonstrated to have great influence on the struc- 
ture of the supercritical, quasi-perpendicular shock. If 
a high ujpe/flce ratio is used, a Buneman instability be- 
tween reflected ions and the inflow electrons plays an 
important r ole in the foot electric field i n electron in- 
ertial scale ( Shimada and Hoshino 2000l ). and on this 
scale the electron pressure term may become as impor- 
tant as the Lorentz term and Hall term analyzed in 
this paper. Simultaneously, two-dimensional PIC sim- 
ulations found that a supercritical quasi-perpendicular 
shock may be emitted waves, such as oblique whistler 
waves and lowe r -hybr id waves within the shock front 



2009). which may also change the 
structures of the shock. In addition, Alfven waves 



( Lembege et al 



( Lu et alJ 120091 ) or upstream turbulence generated by 
summi ng many individual circularly polarized Alfven 
waves (Giacalonc and Deckerll2010l: iGuo and Giacalond 
2010) are always observed in the solar wind. Clear rip- 
ples can be identified at the front of the shock after 
the interaction of upstream turbulence structures with 
a shock. How they will affect the cross shock electric 
field/potential and the mechanisms of ion acceleration 
at perpendicular shocks is our future topic. 
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Fig. 1 The time evolution of B y versus X. Thick curves A, 
B, and C indicate the shock profiles that are chosen in the 
stack plots within one cyclic self-reformation of the shock 
front (with 25% pickup ions). 
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Fig. 2 (a) phase space plots (Vi X — X) of the solar wind 
ions at the shock profile A (t =10.4fi^ ). And the cross 
shock potential is also shown for reference (in blue), (b) 
the density profiles of the pickup ions, solar wind ions and 
their total are in blue (line 1), red (line 2) and green (line 
3), respectively. The main magnetic field B y (black curve) 
is also shown for reference. "O" denotes the position of the 
old ramp (marked by a vertical dashed) during one reform- 
ing cycle, (c) the cross shock electric field E x is indicated 
by the black curve. The additional curves indicate the con- 
tributions to E x of the Hall terms on the right of Eq. ((3}: 
pickup ion Hall term -^-(Jpis x B)^ (line 1), solar wind ion 
Hall term -^-(Jsws X B)^ (line 2), and their total (line 3). 
(d) corresponding plots for the Lorentz terms on the right 
of Eq. (|3|: pickup ion Lorentz term — ^^KUpjs x ~B) X 
(line 1), solar wind ion Lorentz term — ns J Vs (JJsws x H) x 
(line 2), and their total (line 3). The cyan, yellow and ma- 
genta highlight the regions: stationary foot, nonstationary 
foot and the nonstationary ramp, respectively. 
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Fig. 3 (a) phase space plots (Vi X — X) of the solar wind 
ions at the shock profile A (t — 10.9^1^). And the cross 
shock potential is also shown for reference (in blue), (b) the 
density profiles of the pickup ions, solar wind ions and their 
total are in blue (line 1), red (line 2) and green (line 3), 
respectively. The main magnetic field B y (black curve) is 
also shown for reference. "O" and "N" denote the positions 
of the old and new ramps (marked by a vertical dashed 
lines) during one reforming cycle, (c) the cross shock electric 
field E x is indicated by the black curve. The additional 
curves indicate the contributions to E x of the Hall terms on 
the right of Eq. pickup ion Hall term ^-(Jpi s x B)a- 

(line 1), solar wind ion Hall term — (Jsws x B)» (line 2), 
and their total (line 3). (d) corresponding plots for the 
Lorentz terms on the right of Eq. ((3]): pickup ion Lorentz 
term — npls (Upis x TZ)x (line 1), solar wind ion Lorentz 
term - (Usws x B) x (line 2), and their total (line 
3). The cyan, yellow and magenta highlight the regions: 
stationary foot, nonstationary foot and the nonstationary 
ramp, respectively. 




Fig. 4 (a) phase space plots (Vi X — X) of the solar wind 
ions at the shock profile A (t =11.3fl~^). And the cross 
shock potential is also shown for reference (in blue), (b) 
the density profiles of the pickup ions, solar wind ions and 
their total are in blue (line 1), red (line 2) and green (line 
3), respectively. The main magnetic field B y (black curve) 
is also shown for reference. "N" denotes the position of the 
new ramp (marked by a vertical dashed) during one reform- 
ing cycle, (c) the cross shock electric field E x is indicated 
by the black curve. The additional curves indicate the con- 
tributions to E x of the Hall terms on the right of Eq. @: 
pickup ion Hall term ^-(Jp/ S x ~B) X (line 1), solar wind ion 
Hall term -^-(Jsws x B)^ (line 2), and their total (line 3). 
(d) corresponding plots for the Lorentz terms on the right of 
Eq. 0: pickup ion Lorentz term — npla (Upj s x B)^ (line 

1) , solar wind ion Lorentz term — n? ™ s (Usws x B)^ (line 

2) , and their total (line 3). The cyan and magenta highlight 
the regions: stationary foot and the nonstationary ramp, 
respectively. 
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Fig. 5 The time evolution of B y versus X. Thick curve A' 
indicates the shock profile that is chosen in the stack plots 
at the stationary shock with 55% pickup ions. 




Fig. 6 (a) phase space plots (Vi X — X) of the solar wind 
ions at the shock profile A' (t =10.4f2~ 1 ). And the cross 
shock potential is also shown for reference (in blue), (b) 
the density profiles of the pickup ions, solar wind ions and 
their total are in blue (line 1), red (line 2) and green (line 
3), respectively. The main magnetic field B y (black curve) 
is also shown for reference, (c) the cross shock electric field 
E x is indicated by the black curve. The additional curves 
indicate the contributions to E x of the Hall terms on the 
right of Eq. ((3J): pickup ion Hall term -^-(Jpi s x B)^ (line 
1), solar wind ion Hall term ^-(Jsws X B)^ (line 2), and 
their total (line 3). (d) corresponding plots for the Lorentz 
terms on the right of Eq. (J3]): pickup ion Lorentz term 
— npls (Up/s x B)j; (line 1), solar wind ion Lorentz term 
_"Ms^(lj SWs x (i ine 2), and their total (line 3). The 
cyan and magenta highlight the regions: stationary foot and 
the nonstationary ramp, respectively. 



11 



References 

Burgess, D., Wilkinson, W.P., Schwartz, S.J.: J. Geophys. 
Res. 94, 8783 (1989) 

Burgess, D., Lucek, E.A., Scholer, M., Bale, S.D., Balikhin, 
M.A., Balogh, A., Horbury, T.S., Krasnoselskikh, V.V., 
Kucharek, H., Lembege, B., Mobius, E., Schwartz, S.J., 
Thomsen, M.F., Walker, S.N.: Space Sci. Rev. 118, 205 
(2005) 

Chapman, S.C., Lee, R.E., Dendy, R.O.: Space Sci. Rev. 
121, 5 (2005) 

Dunlop, M.W., Balogh, A., Glassmeier, K.-H., Robert, P.: 
J. Geophys. Res. 107, 1384 (2002) 

Eastwood, J.R, Bale, S.D., Mozer, F.S., Hull, A. J.: J. Geo- 
phys. Res. 34, L17104 (2007) 

Giacalone, J., Decker, R.: Astrophys. J. 710, 91 (2010) 

Gloeckler, G., Geiss, J., Roelof, E.C., Fisk, L.A., Ipavich, 
F.M., Ogilvie, K.W., Lanzerotti, L.J., Steiger, R.von, 
Wilken, B.: J. Geophys. Res. 99, 17637 (2002) 

Guo, F., Giacalone, J.: Astrophys. J. 715, 406 (2010) 

Hada, T., Oonishi, M., Lembege, B., Savoini, P.: J. Geo- 
phys. Res. 108, 1233 (2003) 

Heinz, H., Paul, W., Binder, K.: Phys. Rev. E 72, 066704 
(2005) 

Hcllinger, P., Travm'cek, P., Lembege, B., Savoini, P.: Geo- 
phys. Res. Lett. 34, L14109 (2007) 

Kucharek, H., Scholer, M.: J. Geophys. Res. 100, 1745 
(1995) 

Lee, M., Ip, W.-H.: J. Geophys. Res. 92, 11041 (1987) 
Lee, R.E., Chapman, S.C., Dendy, R.O.: Ann. Geophys. 23, 
643 (2005) 

Lembege, B, Dawson, J.M.: Phys. Fluids 30, 1767 (1987) 
Lembege, B., Savoini, P.: Phys. Fluids 94, 8783 (1992) 
Lembege, B., Giacalone, J., Scholer, M., Hada, T., Hoshino, 
M., Krasnoselskikh, K., Kucharek, H., Savoini, P., Tera- 
sawa, T.: J. Geophys. Res. 110, 161 (2004) 
Lembege, Savoini, P., B., Hellinger, P., Travm'cek, P.: J. 

Geophys. Res. 114, A03217 (2009) 
Liewer, P.C., Goldstein, B.E., Omidi, N.: J. Geophys. Res. 

98, 15211 (1993) 
Lipatov, A.S., Zank, G.P.: Phys. Rev. Lett. 82, 3609 (1999) 
Lobzin, V.V., Krasnoselskikh, V.V., Bosqued, J.-M., Ping- 
on, J.-L., Schwartz, S.J., Dunlop, M.W.: Geophys. Res. 
Lett. 34, L05107 (2007) 
Lu, Q.M., Hu, Q., Zank, G.P.: Astrophys. J. 706, 687 
(2009) 

Lu, Q.M., Wang, R.S., Xie, J.L., Huang, H., Lu, S., Wang, 

S.: Chin. Sci. Bull. 56, 1174 (2011) 
Matsukiyo, S., Scholer, M., Burgess, D.: Ann. Geophys. 25, 

283 (2007) 

Matsukiyo, S., Scholer, M.: J. Geophys. Res. 116, A08106 
(2011) 

Mazelle, C, Lembege, B., Morgenthaler, A., Meziane, K., 
Horbury, T.S., Genot, V., Lucek, E.A., Dandouras, I.: 
AIP Conf. Proc. 1216, 471 (2010) 

Mobius, E., Hovestadt, D., Klecker, B., Scholer, M., Gloeck- 
ler, G., Ipavich, F.M.: Nature 318, 426 (1985) 

Mobius, E.: Adv. Space Res. 6, 199 (1986) 

Neugebauer, M., Lazarus, A.J., Balsiger, H., Fuselier, S.A., 
Neubauer, F.M., Rosenbauer, H.: J. Geophys. Res. 94, 
5227 (1989) 



Nishimura, K., Matsumoto, H., Kojima, H.: J. Geophys. 

Res. 108, 1182 (2003) 
Oka, M., Terasawa, T., Noda, H., Mukai, T., Saito, Y.: 

Geophys. Res. Lett. 29, 1688 (2002) 
Omidi, N., Winske, D., Wu, C.S.: ICARUS 66, 165 (1986) 
Omidi, N., Winske, D.: J. Geophys. Res. 92, 13409 (1986) 
Pritchett, PL.: J. Geophys. Res. 106, 3783 (2001) 
Quest, K.B.: Phys. Rev. Lett. 54, 1872 (1985) 
Richardson, J.D., Kasper, J.C., Wang, C, Belcher, J.W., 

Lazarus, A.J.: Nature 454, 63 (2008) 
Sagdeev, R.Z.: Rev. Plasma Phys. 4, 23 (1966) 
Sagdeev, R.Z., Shapiro, V.D., Shevchenko, V.I., Szego, K.: 

Geophys. Res. Lett. 13, 85 (1986) 
Scholer, M., Matsukiyo, S.: Ann. Geophys. 22, 2345 (2004) 
Shapiro, V.D., and Ucer, D.: Planet. Space Sci. 51, 665 

(2003) 

Shimada, N., Hoshino, M.: Astrophys. J. 543, L67 (2000) 
Shimada, N., Hoshino, M.: J. Geophys. Res. 110, A02105 
(2005) 

Shinohara, I., Fujimoto, M., Takaki, R., Inari, T.: IEEE 

Trans. Plasma Sci. 39, 1173 (2011) 
Vasyliunas, V.M., Siscoe, G.L.: J. Geophys. Res. 81, 1247 

(1976) 

Walker, S.N., Alleyne, H.St.C.K., Balikhin, M.A., Andre, 

M., Horbury, T.S.: Ann. Geophys. 22, 2291 (2004) 
Wu, P., Winske, D., Gary, S.P., Schwadron, N.A., Lee, 

M.A.: J. Geophys. Res. 114, A08103 (2009) 
Yang, Z.W., Lu, Q.M., Lembege, B., Wang, S.: J. Geophys. 

Res. 114, A03111 (2009a) 
Yang, Z.W., Lu, Q.M., Wang, S.: Phys. Plasma 16, 124502 

(2009b) 

Yang, Z.W., Lembege, B., Lu, Q.M.: J. Geophys. Res. 116, 
A10202 (2011a) 

Yang, Z.W., Lembege, B., Lu, Q.M.: J. Geophys. Res. 116, 
A08216 (2011b) 

Zank, G.P., Pauls, H.L., Cairns, I.H., Webb, G.M.: J. Geo- 
phys. Res. 101, 457 (1996) 



This manuscript was prepared with the AAS IATr^X macros v5.2. 



